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THE FRICTIONLESS FLOW 111 THE REGION AROUND TWO CIRCLES* 

By M. Lagally 



1 • Obje c t^_,of ...t.he__i nve^,ti^ T^o noninter- 

secting circles are given in a plane; a frictionless flow is 
sought which fills the entire space A around "both of the cir- 
cles* The flow must he free of singularities in finite space; 
however, it may possess a velocity at infinity. Also if this is 
given, the flow is therewith nonuniquoly determined; there re- 
sults a doubly infinite nuaiber of solutions to the problem, the 
first of which will be unique provided that we write down a de- 
termined value for the circulation around each of the two cir- 
cles. The solution of the nroblom is carried out with the aid 
of the conformal representation of the space A as a rectangle; 
this gives four simple flows from which the general flow appears 
to bo linearly built up. Of particular interest is the determi- 
nation of the flux that passes through the gap between the two 
circles, with which the flow coming from infinity and its de- 
pendence on both of the circulations and the velocity at infini- 
ty are interrelated. 

The investigation is ca.rricd thus far« It attempts to sur- 
pass the boundaries of pure mathematical interest in that it 
possesses as an example, a flow investigation in a multiply- 
connected region. Then the results appear to be carried out by 
means of an appropriate conformal representation of the region 
around two chosen closed curves. Thus we have the basis of an 
exact plane theory of the biplane. Certainly, the known diffi- 
culties of the accomplishment of this representation depend upon 
the region around the given curves* It is easier to choose the 
transformation function so that wo obtain the transformed curves 
which, as contours of the wings of a biplane are at least useful 
to a certain extent, and correspond to 'the Joulcowski monoplane 
wing. Both wings have sharp trailing edges and, moreover, in 
order to obtain a smooth flow at the trailing edges, we have to 
locate a branch point of the transformation function which falls 
at a stagnation point of the flow around the circles.. This ap- 
plication is only pointed out here and will not be entered into 



*"Die reibungslose Stromung im Aussengebiet zweier Kreise." 
Zeitschrift fur angewandte Matheraatik und Mechanik, Vol, 9, No. 
4, August 1929, pp. 299-305. 



2 



IT.A.C.A, Technical Memorandum IIq. 626 



in what follows.* 

2 • Bipolar coordinat e s_and_as so c 

tion.- "ITe make two points (c,o) and ^2 (-c,o) on the x- 

axi s in the z^-plane, the origins of trro polar coordinate systems 

r^, and r^ , so that we may take the complex variable 

z = X + iy for a point P (x,y) in "both the forms: 

z = c + r^^ e^^^ and z = - c + r2 e^ . (l) 

Then the quotient 

= £3 ei - ^) ■ (2) 

z - c ri 

r^ 

contains in a known manner, the "bipolar coordinates and 
1^2 - -i^i . (?ig. la.) ^ 

The curves — = constant are according to the theorem of 

Apollonius, the circles of a circular pencil which contains the 
points and as real zero circles. The curves ^ 

constant are, "by means of the theorem on the peripheral angle, 
circles of a circular pencil v/hich have the points and 
as real points of intersection; each of these circles is split 
into two arcs "by and Q2 for which the corresponding pe- 

ripheral angles 1^2 - ^^1 and - ''^2 total 2Tr.. Both the cir- 

cular pencils are orthogonal. 

^o introduce a now complex variable Z = X -f i Y "by means 
of the equation 

Z = logof-±-°- (3) 

which separates into 

X - loge~, Y = ^2 - ^1 (4) 

*Note "by the Editor: LIr . Stephan Bergmann has called to my at- 
tention, that Julius Bonder has handled the same problem in his 
thesis (Proceedings and reports of the Warsaw Poly t echni cal Soci- 
ety, Vol. 17, 1, 1925). J. Bonder applies the method of images, 
however, he does not recognize that the series found "by him arc 
only partial fractional developments of elliptic functions. (See 
the report of St. Bergmann on "l^ouoro polnischo Arbeiten ubcr Me- 
chanik," this journal. Vol. IX, 1929, pp. 241-247.) furthermore, 
Dmitri Ri ab '^uchinsky has treated in his thesis (1922) in part 2, 
Chap.. IV ( "Llouvement cy clique ;plan d»un liquide autour d*un s el- 
ide qui^ se liiout parallelement a une parol rectiligne) a somewhat 
special problem with the aid of Jacobi functions. 
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so that ue ottain a conforraal representation of the orthogon^vl 
circular pencils. (Fig. Ih.) 

The circles ~ = constant correspond in a unique raanner 

to the straight lines X = constant; in particular, the zero 
points and Qi to the lines X = - oo and X = + cc. The 

circular arcs 1^2 - ^ constant "between Q, 1 and Q2 corre- 
spond to the lines Y = constant a.nd indeed since the angles 

and ^2 are determined only up to multiple values of 2 n, 
each arc ^2 - -^i = constant, corresponds to infinitely many 
parallels Y. = constant of difference 2 tt. 

The whole z-plane will thus have infinitely many values on 
the Z-plane, "but a strip of width 2 tt hounded "by two paral- 
lels to the X-axis is conformally represented in a unique manner. 
For the following investigation, the strip of the Z-plane lying 
between the lines Y = - rr and Y + tt will he considered as 
the representation of the z-plane cut along the distance Q^. 
By this transformation the point z = co corresponds to Z 0. 

3 . Jhe flow-field a nd its confo rmal reioresent ati on . - (Fig. 
2a, h.)- It will now he necessary to investigate the flow in the 
z-plane which occupies the region A enclosing two circles Ki 
and K2 of the circular pencil with the real zero circles Oi 
and Q2. Let Ki and K2 he doterminod hy means of £2. = \ 

rp ^1 
and ~ = [i,, respectively; thus, if the zero circle lies in 

the interior of and the zero circle Q2 in the interior of 

K2 , then X > 1, |i < 1 . If to shorten, we put log^ \ = a, 
logg |uL = - P , then the circles Ki and K2 correspond to the con- 
formal representation of the two lines X a and X = - P. 

The representation of the region A is thus a rectangle 
which is formed hy the four lines X - a, X = - p, Y = tt, 
y = « TT, The two sides Y-+tt, Y--tt correspond to the 
lower and upper edges, respectively, of a cut along the x-axis 
which is drawn hetv/een the two circles. The rectangle contains 
the point Z - 0 as an inner point. 

If the flow in the z-plane has the circulations ^1 and 
V2 ahout the circles Ki and K2, it. thus keeps as its com- 
plex potential at infinity: 

Q = - i — logez. (5) 

Fx + is the circulation around every closed curve surround- 

ing Ki and K2; the infinitely distant point of the flow it- 
self has the character of a point vortex of circulation 
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r = - (Pi + Ts) . (6) 

By means of the conformal representation, there corresponds to 
it a point vortex with the same circulation at the point Z = 0. 
It accordingly retains as its complex potential in the neighhor- 
hood of this point 

n = - i 2^^- loggZ . (7) 

This appears to "be confirmed "by means of calculation. T7e sul)- 
stitute in (5) "by means of (3) 

z = c 4--^' 
e^ - 1 

thus vre ootain a complex potential 




The function -^--±-1 has a pole of the first order at the point 
e^ - 1 

Z = 0 ; i t i s 

. = 2 + I + ,3 + a3 + ... = 

e2 - 1 Z 6 

; J [2 + |- + as + as Z"" + . . . .] • (8) 

"by the omission of a regular flow in the neighhorhood of Z = 0, 
we ohtain the equa^tion (7) , 

The circulations Ti and" around "both the circles Ki 

and K2 are transf ormed " "by * means of the conformal representa- 
tion into flows Ti and along the sides of the rectangle 
X = a and X = - p, which correspond to the two circles. If 

and are positive, then they occur in the negative direc- 

tion around the rectangle. 

The flow in the 2-plane must "be free of singularities in the 
finite portion of A; however, a velocity at infinity Wod ~ Uco 
- iVoo may he permitted. The resulting translatory flow is ob- 
tained as . 

Q = Wcx> z (9) 

thus the transformed flow in the neighborhood of the point 
Z 0 as 
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Q = c (10) 
- 1 



or after the omission of a regular flow by means of (8) as 

n-^^- . (10«) 

It possesses a doublet at the point Z = 0, 

Thus we can set up the complex potential Q and the com- 
plex velocity ^ ~ ^-^ transformed flow, in the neighbor- 
hood of the point Z = 0, If P (Z) represents a function, 
which is regular in the neighborhood of the zero point, then by 
means of (?) and (lO»)f it is 

p 2 W-„ c . . 

Q = - i logeZ + ~2-— + P (z) ; 

p 2 Wcr. C .V XV 

W - - i ~i;-77 + z)... 11) 

• 2 TT Z Z 

At all other points of the transformed rectangle Q and W are 
regular . 

4* Formation of the complex potential.- The complex poten- 
tial of the transformed flow is determined according to (ll) ex- 
cept for a function which is regula,r in the transformed rectan- 
gle. The complete determination succeeds in the following man- 
ner. The complex velocity in the z-plane is 

^^dfi^dildZ^^^ 1 - 1 ] 

dz dZ dz z + c z - c * 



Since w is a single-valued function in A, the same applies 
to W; that is, W must as a function of Z possess the peri- 
od 2 TT i , thus it has the sajne value at all the transformed 
points of z. Further, since Z = a and X = - p are stream- 
lines, therefore all the same singularities of W situated with- 
in the transformed rectangle and beyond must be symmetrically 
located with respect to both these lines. Thus the whole flow 
in the transformed rectangle will be repeatedly rofloctod over 
both these linos. In every case W must thus possess the peri- 
od 2 (a + p) . 

W is thus a doubly periodic function with the periods 



2 0^1 ^ 2 (a + p) ; 2 oog = 2 tt i 



(12) 
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The periodic rectangle has do^ililG the width of the transformed 
rectangle; for the following considerations a periodic rectangle 
is "basically assumed which is composed of the transformed rec- 
tangle and its reflection on the line X ~ - ^. In it W has 
"besides the singularity at the zero point an additional re- 
flected singularity at the point - 2p. (Pig. 3.) Thus W is 
an elliptic function which may he put in the form 

r 1 . r 1 2 c(uo. - i Voo) 



W = - i i + 



1 



2tt2 2nZ + 2p 

^- + H (Z) (13) 



(Z + 2 p)^ 

w?iere 3 (Z) represents a function whicli is regular in the 
periodic rectangle. This elliptic function i^* determined by 
means of its polos except for an additive constant; we may 
put for hrcvity 

"^03 + i Vcx> = "vvoo with UcD - i Vo^ ~ woot 

so that* 

W = - i [^(Z) - ^(Z + 2p)] - 2c[w^p(Z) - w,^p(Z + 2p)] + i K. 

«o XT 



(14) 



The additive constant is introduced in the imaginary form i K 
for a later evident useful basis; it could be determ.ined by 
means of the substitution of the value of the velocity at a 
point of the flow field; K is thereby proven to be real; 
however, a better way for the determination of K will be 
shortly pointed out. 

The complex potential will be 

Q = - i loge ^' ^^^^ ^^^^ " t(Z+2!3)] + iKZ + K' 

(15) 



Here k» is a new additive constant which can be chosen as we 
please. For simplicity we shall put '^^ - 0. However, it 
shou.ld be noted tliat ^ is also an infinitely mul t ipl e- valued 

function; if logp — changes by ± 2 'rr i by a revolu- 

' ® a (Z 4- 2p) 

tion around a pole, then Q itself changes by ± i. 



*Burkhardt-?aber , "Elliptische Funhtionen," section 30. 
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In order to make Q single-valued in the periodic rectan- 
gle, we must prevent the revolution arou2^d both the poles 0 
and -2 P; this is done "by the introduction of a "branch-cut 
along the combined segments of the x-axis. 

5 • Relati on ship between circulation and flux , - The fun ca- 
tion Q(Z) is an elliptical transcendental one and changes by 
an ^additive constant if the argument of Z increases by a peri- 
od 2 or 2 oOg, written combined as 2 O). 

In order to calculate these changes, Cl(Z) must be put in 
the following form, in which it is to be noted that the factor 
of UcD is an elliptic function: 



Q(Z) = - i -C^ loge---?-in__^ + 2 c uco [C(Z) - ^(Z + 2 p) f 
- 2 i c [^(Z) + t(Z + 3 p)] + i K Z . 



(16) 



Then 



[2(2+20)) -U.Z) - i lose a(Z+2p+203)a(z) 

- icvc» [t(Z+2'a)) + ^(Z+2p+2a))-^(Z)-^(Z+2p)] +2 i k o) . 
From the properties of the t, and o fvinctions* 

az + 2 r^) - C(Z) = 2 7., aiZ^±-.fl^-i ^ 3-2T>.(Z+a)) 

a(z) 

where rj is written for "H, = ^ ( ^i) j ^2 ~ ^ (^'^^) combined, it 
follows immediately for both the periods written separately 

Q(Z -f 2 coj « Q(Z) ^= i [S/tt T g ri^ - 8 c Vo3 'P.i + 2i^a)i] , (17a) 
C(Z -f 2 0:2) - Q(Z) - i [2/n T p ri^ - 8 c Voo r\^^ + 2^0)2] . (i7b) 



It should be noted here that oOi and T]i are real, ^2 and 
r\p^ purely imaginary, that therefore the right side of (17a) is 
purely imaginary, that of (17b) is real^ This remark brings in 
the hydrodynamical meaning of the result. We decompose the com- 
plex potential Q = cp + i V into real and imaginary parts; 
thus cp represents the velocity potential, ^|^ the stream func- 
tion. We connect the points Pi and P2 in a region in which 
0 is single-valued, by any path we please, and decompose the 
difference Q2 - Qi in the same way: 

*Burkhardt-?aber , sections 25-25, 
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so that CP2 - cpi is the flo?7 ^long the path PiP,, and 

^2 - ^1 is the flux through this path, calculated .as positive 

from the right , 

Thus Q(Z + 2 oji) - Q(z) is, except for the factor i, 
the flux throu^;h t'iie horizontal sides ox the periodic rectangle. 
We put n(Z + 2 - C(Z) = ^ 2 i so that ? is the flux 

through the horizontal sides of the transforned rectangle, cal- 
culated as positive from the left, or also, "because of the in- 
variance of Q under a conforiao,l transformation, the flux 
through a chosen line joining two points of the circles and 
Kg, calculated as positive from the loft. 

Similarly C:(Z + 2 '^2) - C^\Z) is the flow along one 6f 
tho'vortical sides of the periodic rectan^'lo, hence - P^; or. 
except for the sign, the circVi.la.tion around the circle Zi. 
(?ig. 3.) 

Thus it follows from (17a, h) 

- 2 i ? ^ i [2/17 r p T;i - 8 c Vc^ 'P.i 'h 2 K ooi] , (18a) 

^" = i [2/tt r P r,2 8 c Vo. + 2 « ^2] • (18b) 

Each of these equations rould permit the calculation of the con- 
stant with the aid of the q\iantities Ti, F, F, Voo , on each 
of which the flow around the two circles depends. We eliminate 
K on the other side of hoth of the equations "by the application 
of the Legendre relation* 

'Hi - '^.2 ^1 = -"2^ , (19) 

60 that we obtain an interesting relation between the hydrody- 
namical quantities : 

2 i P 0)2 + Pi 0)1 = - r P -I- 4 rr c Vro. 

This is simplified by means of (6) and (12) into 

2 TT'P - 4 rr c v^ - a + (20) 

This equation gives the flux passing between the two circles if 
the velocity at infinity and the circulations around both the 
circles are given. 



*Burkhardt-Paber , section 25, (14). 
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S« Si^-QIlr^-Si on_of _the_2oj^ihl^ If the opposite con- 

ditions for two circles are g:iven, hy means of the three con- 
stants a, p, and c, then the flow in the region surrounding 
"both the circles is therewith nonuniquely determined; likewise, 
also if we have given the velocity at infinity w^^ = Uq^ - i r^. 
According to (16) a, doubly infinite number of different flows 
occur by means of the choice of the constants P and 
also, if we take from (18b), 

by the choice of the constants F and r^. 

Wo introduce further, besides P ~ - {T F^) , the quanti- 
fy 

i = 1 1 ^ l2 

thus putting 

i i / 4 i 3 ri^N a 

i K = ~ r - ~ i 1 r + ^ c Voo Ti^ . 

4tT 4tT V TT / ^ IT 00 12 ^ 

so that the complex potential (IG) appears in the following 
form (out of which the unsymmotric advantages of the quantity 
P before a could be, set aside by means of (12) and (19) ): 

+ 2 c Ucc Caz)- az+2p) ] -2icv,, [az)+ r.Z+2p)+ --^^^ Z ] . j 

The four-term aggregate on the right side a.ppears to show that 
the whole flow arises from the linear supposition of four sepa- 
rate flows, for each of which onl^^ one of the four quantities 
r t r, Uoc, and Voo, is different from zero. 

Bv the first partial flow -i- Z, both the circulations 

are equal and- cj^jposi te , if we return to the z -plane: = - 

= V /2. ThG streamlines are the circles of the pencil; the 
velocity at infinity is zoro« 

By the second partial flow with the factor P, the two 
circulations are equal: p^ P^ = - p/2; the velocity at infin- 
ity is zero * 

The third partial flow with the factor u^^, is a noncircu- 
latory translatory flow in the direction of the line of centers 
of both the circles. 
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The fourth partial flow with the factor is a non- 

circulatory tr?.nslatory 'flow perpendicular to the direction of 
the line of centers of hoth the circles. 

If we introduce the constants • and in (20), then 

we may recognize in this equation the portion of the partial 
flows in the flux passing through the gap hetwoen the two circles. 

2ttP - 4ttcv^ ^-:^^r - ^^^r' . (22) 



For a noncirculatory flow, (22) redu.ces to 



F = 2 c V. 



(22') 



this is also the flux of an unint er rujjt ed flow through a canal 
of width 2c; the flux v/hich passes "between the two circles is 
independent of the actual distance of both the circles, and 
does not increase if we replace "both of the circles hy others 
of the sane pencil; the effective v;idth of the gap "between two 
circles is the distance of the zero circles of the pencil de- 
terrained hy the two circles. 



• l*lj-.^Q-^*-'^<^!^-lQ^ Qf Z^""^^. radii and c ent or- t o- c ent er distance 
2.^._-^2.t'^— • ^'^"^^ w^o of "bipolar coordinates has the disad- 
vantage that "both the circles and are not given hy 



their radii and and the abscissas a 



and -a. 



of 



a^ + ag > 



1 '^^"^ ""2 

their centers, or the cent er- t o-cent er distance d 

hut separately hy three constants a, P, and c, which de- 
termine the circles only indirectly. It remains only to add 
the equations which serve to introdi\ce R^, Ro and d. 



(Fig. 



^a 



• ) 



The equation of the first circle 



K ^ in bipolar co- 



ordinates is 



^2 



.a. 



ordinates we can replace 



in order to bring in rectangular co- 
Iz + c! 



r o 

~ by 
^1 



Ci 



A simple trans- 



formation give's the equation of the circle 



- 2 X c coth (X + 



0, 



from which 



= c coth a; 



R 



____c 

s i nh a 



appears to be taken. 
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For the second circle, , there follo77S from — ^ = e » 

^ 1 

the corresponding valiies 

- = - c coth j3; Ho ^ ~t-^ — o . 

^ Sinn p 

Thus if "both the radii R^^,., and the center-t o-contor 

distance d are given, we have the three equations 

= —r^r- — ; R2 = — rl = c(coth a H- coth B) (23) 

^ smh a , smh p ^ 

for the calculation of a, p, and c. From these three simple 
and admirably detached equations, the introduction of any (more- 
over conformally invariant) auxiliary quantity is permitted, 
namely, the intersecting;; angle of "both the circles. 

In general the intersecting; an/<;le of two circles is given 

- R " - 

COS (jj = ; 

2 R.E^ 

for d >R^ + Eo, the an^^le will "be imaginary. We put a) = i J 
so that cos u) = cosh J, thus we have for the calculation of 
the auxiliary quantity J, the equation 

d — E - K 

cosh J = ~ ~ . (24) 

2 H:,E2 

J appears also to "be expressed "by the system of quantities a, 
p, and c; from (23) and (24) it follows that 

cosh J = cosh (a + p); thus J = a + p. 

Now this gives the solution of the equation (23) in simplest 
form 

EE E E 

c = ~— - sinh J; sinh a ^ — sinh J, rinh P = ~ sinh J (25) 
d d d 

where J is to "be taken from (24) • 



Translation "by liiltcn J. Thompson, 
University of Michigan, Ann Arbor, Michigan. 
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Figs. 2a, 213,3 





